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ABSTRACT 

If core collapse leads to the formation of a rapidly rotating bar-unstable proto-neutron star surrounded 
by fall-back material, then we might expect it to cool and fragment to form a double (proto)-neutron 
star binary into a super-close orbit. The lighter star should survive for awhile, until tidal mass loss 
propels it toward the minimum stable mass of a (proto)-neutron star, whereupon it explodes. Imshennik 
& Popov have shown that the explosion of the unstable, cold star can result in a large recoil velocity of 
the remaining neutron star. Here, we consider several factors that mitigate the effect and broaden the 
range of final recoil speeds, in particular the finite velocity and gravitational deflection of the ejecta, a 
range of original masses for the low mass companion and its cooling history, rotational phase averaging 
of the momentum impulse from non-instantaneous mass loss, and the possibility of a common envelope 
phase. In spite of these mitigating factors, we argue that this mechanism can still lead to substantial 
neutron star recoil speeds, close to, or even above, 1000 km s~^. 

Subject headings: stars: neutron - supernovae: general - pulsars: general - stars: rotation 



1. INTRODUCTION 

It is recognized that radio pulsars have peculiar space 
• velocities between « 30kms~^ and « 1600 km s""'^, signif- 
' icantly greater than those of their progenitor stars. The 
. highest speed ever recorded is from the Guitar Nebula pul- 
sar (Cordes et al. 1993; Cordes & Chernoff 1998), while 
the lowest has been recently measured for B2016-I-28 by 
Brisken et al. (2002) from very accurate VLBA pulsar par- 
allaxes. Statistical studies, aimed at inferring the peculiar 
, velocity at birth from the observed speed, give mean three- 
' dimensional velocities of 100 - 500 kms~^ for the isolated 
\ pulsars (Lyne & Lorimer 1994; Lorimer, Bailes, & Har- 
. rison 1997; Hansen & Phinney 1997; Cordes & Chernoff 
' 1998; Arzoumanian, Chernoff, & Cordes 2002). 

An early explanation for the large space velocities called 
for recoil in a close binary that becomes unbound at the 
'■ time of (symmetric) supernova explosion (Blaauw 1961; 
. Iben & Tutukov 1996). Now, a number of observations 
' hint at a natal origin of these high space velocities, or at 
a combination of orbital disruption (when the progenitor 
lives in a binary) and internal kick reaction. Evidence for 
a kick at the time of neutron star birth is now found in a 
variety of systems: in runaway 0/B associations (Leonard 
& Dewey 1992), in highly eccentric Be/NS binaries (van 
den Heuvel & Rappaport 1986; Portgies-Zwart & Verbunt 
1996), in the binary pulsar J0045-7319 (Kaspi et al. 1996; 
to explain its current spin-orbit configuration) , and in dou- 
ble neutron star binaries such as B1913-I-16 (Bailes 1988; 
Weisberg, Romani, & Taylor 1989; Cordes, Wasserman, & 
Blaskiewicz 1990; Kramer 1998; Wex, Kalogera & Kramer 
2000) where misalignment between the spin and orbital 
angular momentum axes indicates velocity asymmetry in 
the last supernova. All these observations support the 
view that the formation of neutron stars is accompanied 
by anisotropic explosion. This notion is bolstered by evo- 
lutionary studies of binary populations (Dewey & Cordes 



1987; Fryer & Kalogera 1997; Fryer, Burrows, & Benz 
1998) and by studies (e.g. Cordes & Wasserman 1984) on 
the survival of binary systems into their late evolutionary 
stages after supernova explosion. 

Among the physical processes that have been proposed 
to account for the kicks are large-scale density asymmetries 
seeded in the pre-supernova core (leading to anisotropic 
shock propagation), asymmetric neutrino emission in pres- 
ence of ultra-strong magnetic fields (see Lai, Chernoff, & 
Cordes 2001 for a review), or off-centered electromagnetic 
dipole emission from the young pulsar (Harrison & Tade- 
maru 1975). However, none of these mechanisms can ex- 
plain kicks as large as ~ 1600 km s~^. 

In this paper we reconsider the idea first put forward by 
Imshennik & Popov (1998) that, in the collapse of a rotat- 
ing core, one or more self-gravitating lumps of neutronized 
matter may form in close orbit around the central nascent 
neutron star, transfer mass in the short lived binary, and 
ultimately explode, causing the remaining, massive neu- 
tron star to acquire a substantial kick velocity, as high 
as the highest observed. The light member explodes as 
mass transfer drives it below the minimum stable mass for 
a neutron stars. In the light star, stability is lost upon 
decompression by the /3— decaying neutrons and nuclear 
fissions by radio-active neutron-rich nuclei (Colpi, Shapiro 
& Teukolsky 1989, 1991; Blinnikov et al. 1990), that de- 
posit energy driving matter into rapid expansion (Colpi, 
Shapiro & Teukolsky 1993; Sumiyoshi et al. 1998). The 
kick tries its origin from the orbital motion of this evanes- 
cent super-close binary, that forms in the collapse of a 
rapidly rotating (isolated) iron core. 

We will study several effects that may modify the mag- 
nitude of the kick, such as gravitational bending of the 
exploding debris, rotational averaging of the momentum 
impulse, orbit decay, and delayed neutron star cooling. 

Formation of a proto-neutron star companion around 
the main neutron star has never been verified in numerical 
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simulations due to computational limitations. For this rea- 
son wo elaborate on a study of Bonnell (1994) on the for- 
mation of binary /multiple systems in collapsing gas cloud 
cores and its extension to the stellar core collapse in the 
aftermath of a supernova explosion (BonncU & Pringle 
1995) to motivate our working hypothesis. The forma- 
tion of such exotic binaries has been conjectured to occur 
in many works (Ruffini & Wheeler 1971; Clark & Eard- 
ley 1977; Blinnikov, Novikov, Perevodchikova, & Polnarev 
1984; Nakamura & Fukugita 1989; Stella & Treves 1987). 



2. LIGHT FRAGMENTS AROUND PROTO-NEUTRON STARS 

2.1. The Scenario 

Formation of a light companion around a main body 
implies breaking of spherical and axial symmetry dur- 
ing collapse and following core bounce. During dynam- 
ical collapse, unstable bar modes (m = 2) can grow in a 
fluid (even non rotating) that may end with fragmenta- 
tion. However this is known to occur only if the cloud 
core contracts almost isothermally, as in the case of star's 
formation from unstable cold gas clouds (Bonnell 1994). 
Core collapse in type II supernovae is far from isothermal 
(it is described by a effective polytropic index 7 ~ 1.3) 
so that instabilities of this type do not have time to grow 
(Lai 2000) ^, and simulations of non-axisymmetric rotat- 
ing core collapse confirm this trend (Rampp, MuUer & 
Ruffer 1998; Centrella, New, Lowe, & Brown, 2001). Can 
fragmentation/fission be excited after core bounce ? 

Rapid rotation in equilibrium bodies is known to ex- 
cite non-axisymmetric dynamical instabilities and these 
instabilities may grow in the proto-ncutron star core. In- 
terestingly, core collapse simulations of unstable rotating 
iron cores (Heger, Langer & Woosley 2000; Fryer & Heger 

2000) or polytropes (Zwerger & Muller 1997) indicate that 
proto-neutron stars, soon after formation, can rotate dif- 
ferentially above the dynamical stability limit set when the 
rotational to gravitational potential energy ratio Trot/lW^I 
is larger than the value pdyn = 0.25 — 0.26 (Saijo et al. 

2001) . Strong non-linear growth of the dominant bar- 
like deformation (m = 2) is seen in these cores (described 
as polytropes by Rampp, Muller, & Ruffert 1998). How- 
ever there is no sign of fission into separate condensations. 
The bar evolves, producing two spiral arms that transport 
the core's excess angular momentum outwards (see also 
Shibata, Baumgarte & Shapiro 2000). This reduces the 
bar's angular momentum such that a single central body 
is formed. How can a body develop local condensations 
when bar-unstable ? 

According to Bonnell's picture, the evolution of the bar 
instability is more complex, in reality. If the rapidly spin- 
ning proto-neutron star core goes bar unstable when sur- 
rounded by a fall-back disk, then matter present in the 
bar-driven spiral arms interacts with this material. The 
sweeping of a spiral arm into fall-back gas can gather suf- 
ficient matter to condense into a fragment of neutronized 
matter. This occurs because the m = 1 mode grows dur- 
ing the development of the m = 2 mode. The m = 1 mode 
causes the displacement of the unpinned (free-to-move) 
core and creates an off-center spiral arm that sweeps up 



more material on one side than the other during " continu- 
ing" accretion. The condensation eventually collapses into 
a low mass neutron star. 

Detailed simulations confirming or dismissing the occur- 
rence of such instability are still lacking, so further con- 
siderations of the lump masses, temperatures and entropy 
contents are necessarily speculative. 

2.2. Cooling scenarios and the minimum mass 

Here we wish to explore the possibility that a light 
(proto-)neutron star forms around the main central body, 
lives for a while, and later explodes, imprinting a kick to 
the neutron star that remains (due to linear momentum 
conservation) before gravitational waves or hydrodynami- 
cal effects can drive it toward the central star. Can a light 
(proto)-neutron star form from the condensation of ma- 
terial accumulated in the off-centered spiral arm ? What 
limits can be imposed on its mass ? 

Cooling plays a key role in addressing these questions. 
Goussard, Haensel & Zdunik (1998), Strobel, Schaab & 
Weigel (1999), and Strobel & Weigel (2001) have shown 
that the value of the minimum stable mass mmmc for a neu- 
tron star (located at the turning point in the mass-radius 
relation of equilibria) is a function of the temperature: it 
varies from >1Mq at 50-100 milliseconds after core bounce 
(setting the actual value of the mass of the central neutron 
star), to ^ O.7M0 after ^ 1 second, down to ^ 0.3Mq af- 
ter 30 seconds, reaching the value of ~ 0.0925AfQ (Baym, 
Pethick & Sutherland 1971) known for cold catalyzed mat- 
ter (T < 1 MeV) several hundred seconds later. Thus, the 
lump of nuclear matter gathered in the spiral arm by the 
instability may become self-bound if its mass m is above 
the minimum corresponding to that particular tempera- 
ture. Once formed it is stabilized against expansion by 
cooling. The actual value of m is thus determined by the 
overall dynamics of collapse after core bounce and varies 
between 0.0925 and ^ IMq : It depends on the time at 
which the instability sets in, on the amount of fall-back 
material (potential reservoir of matter in the lump) and 
on the cooling history. 

The value of m and of the mass ratio q = m/M in the 
binary (with M the heavier of the two stars) remains un- 
predictable to us at this level, so we may just depict three 
possible scenarios for the formation and evolution of this 
evanescent binary: case (A) when the instability sets after 
few hundreds of milliseconds or a second after core bounce 
and binary formation occurs so that m > O.7M0 (imply- 
ing a pre-existing iron core of large mass if the primary is 
as massive as I.4M0); case (B) when a lighter star forms 
around a main body with m ~ 0.2 — O.4M0 after several 
tens of seconds from core bounce; case (C) when cooling is 
sufficiently advanced that the minimum mass approaches 
its asymptotic value and the binary can have m < 0.2Mq. 

The magnitude of a natal kick in (A) is difficult to es- 
timate and we must await for realistic simulations of core 
collapse. Values oi m/M above a given threshold in a bi- 
nary are known to lead to unstable mass transfer and thus 
to final coalescence. We note that the evolution in this case 
might be similar to that described in coalescing neutron 
star binaries with slightly unequal masses (see Rosswog et 

^Large-scale asymmetries imprinted in the iron core prior to collapse may lead to anisotropic explosions that produce kicks as indicated by 
Goldreich, Lai & Sahrling (1996). Whether they lead to fragmentation is unknown. 
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al. 2000) where it has been shown that large kicks can be 
acquired as a consequence of mass loss via a wind. Case 
(B) and (C) will be explored in this paper. Imshennik & 
Popov (1998) estimated neutron star recoil speeds of order 
1500- 2000 km s^i for scenario (C). 

In cases (B) and (C), one factor that could lower the 
final neutron star speed is the finite velocity of the ejecta. 
Colpi, Shapiro & Teukolsky (1993) have shown that in 
the dynamical phase of the explosion the ejecta can attain 
speeds varying from 10,000 to ~ 50, 000 km s~^ (the upper 
bound being related to the value of the binding energy of 
the star at the minimum mass relative to a dispersed state 
of iron) . These speeds are close to the escape velocity from 
the binary and therefore the final kick imparted to the re- 
maining neutron star may be influenced substantially by 
gravitational deflection of the ejecta, and velocity phase 
averaging during the explosion. A further threat would 
be a rapid decay of the orbital separation due to unsta- 
ble mass transfer that may lead to final coalescence, and 
emission of gravitational waves. 

In §3 we will describe the gravitational bending of the 
ejecta, while in §4 we will study initial conditions and, 
subsequently, orbit decay and mass exchange, for cases 
(B) and (C). We will then give an estimate of the kick 
speed including phase- velocity averaging. 



3. GRAVITATIONAL BENDING OF THE EXPLODING 
DEBRIS AND THE FINAL KICK 

Consider now the instant at which the light secondary of 
mass m explodes, having reached the dynamical instabil- 
ity point (m ~ m,iimc), while orbiting the heavier, primary 
neutron star of mass M. 

Obtaining the final kick speed imparted to the heavy 
neutron star is complicated by several factors. First, some 
of the material ejected from the exploding star can re- 
main bound to the system and eventually accretes onto 
the remaining star. Second, material that escapes may 
not be moving at extremely large speeds, even right at 
ejection, and therefore emerges with a different momen- 
tum than it has at the point of explosion. Third, the mass 
of ejected material may be of order 10% of the total mass 
of the system, and self gravity could play a role. Fourth, 
the ejecta are a fluid, and, depending on how much time 
elapses between neutron star formation and the explosion 
of the low mass star, move through ambient gas left over 
from the original supernova explosion. If the ejecta are 
slowed significantly by the surrounding gas, some material 
that would be judged unbound upon ejection might actu- 
ally end up falling back toward the neutron star, possibly 
to accrete onto it. 

Suppose that relative to the remaining heavy neutron 
star, the orbital velocity of the exploding star is V at the 
instant of explosion (i = 0), and the position of the center 
of mass of the exploding star is at r. Let us write the final 
velocity of the remaining neutron star as 



mV 



(M 



(1) 



where F = If the ejecta are expelled isotropically rel- 
ative to r, and at speeds large compared with V, then we 
expect 77 = |rj| ~ 1, in which case the remaining neutron 



star recoils with the maximum possible kick speed 

mV 

kick.mE 



(m + M) 



(2) 



as implied by momentum conservation in the center of 
mass of the binary. 

To get a more realistic, but still approximate, idea of the 
size of 77 = Vkick /Vkick. max, wc work to zeroth order in the 
ejected mass, and also treat the trajectories of the ejecta 
ballistically; this is equivalent to ignoring the third and 
fourth complications listed above entirely. By working to 
zeroth order in the mass of the ejecta, we can also assume 
that all particles are ejected virtually from a single point, 
r. Let the velocity of a particle relative to the exploding 
star be w, and assume that w is distributed isotropically 
about the position of the exploding star with probability 
P{w). Since we work to zeroth order in the mass of the 
ejecta, we may take the heavy neutron star to remain at 
fixed position in the calculation. 

The problem that remains is the orbital mechanics for 
each particle relative to the massive star. The orbits are 
characterized by the conserved quantities 



E 



1, 



\V + w\' - 



GM 



J = rx{V + w) 

GMr , 
A = h(F + i(j)xJ, 



(3) 



which are the orbital energy and angular momentum per 
unit mass, and the Runge-Lcnz vector, respectively. We 
note that 



A = \A\= GMe = yJ{GMf + 2£;j2; 



(4) 



e < 1 for bound orbits, e > 1 for unbound orbits. 

To calculate 77, we need to find the momentum per unit 
mass carried away by the ejecta that escape to infinite 
distance. First, we need to identify the initial velocities of 
particles that escape. For these, we must have E > 0, and 
from equation (^ we find the condition 



V-w = Vwfi > 



GM (V^ + w^) 



(5) 



where V = \V\ and w = \w\, and /i = V-w/Vw. Sec- 
ond, we need to find the mapping between initial and final 
velocity. To this aim, it is useful to define a coordinate 
system in which 63 = J, ei = A, and 62 = JxA, where 
J ~ J/ J and A = A/A. The orbit is then confined to 
the 1 — 2 plane, and we can read off the final velocity in 
this system from equation (15.12) in Landau and Lifshitz, 
Mechanics: 



Voc ^ [-AV2E + {2E/GM)JxA]. 



(6) 



We wish to express the final velocity iu a coordinate 
system fixed in the binary at the point of explosion. We 
thus choose, explicitly 



r = rx 



and 



V = Vy 



(7) 



so that the orbital angular momentum of the light star 
relative to M just prior to explosion points along z. We 
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have explicitly selected the point of explosion to occur at 
either peri- or apocenter of the orbit, and for simplicity, we 
specialize to circular orbits below. (Both tidal effects and 
gravitational radiation will tend to circularize the orbits, 

but, even for eccentric orbits, we expect tidal disruption 
to be likeliest at pcriccnter.) In this coordinate system, 

J = -yrwz + zr{V + Wy) (8) 

A = x[-GM + r{V + Wyf + rwl] 
+y[-rWa;{V + Wy)] + z{-rw^W;,) (9) 

and 

^_ GM ^ wl + jV + Wy^ + w^, ^ ^^^^ 

The velocity of the ejecta can thus be expressed in terms 
of the known explosion parameters {w, V, and r) by com- 
bining equation (6) with equations (7)-(10). The result 
is that the total momentum carried off per unit mass of 
ejecta is 

u = [ d^wP{w)v^{w), (11) 

JE{W)><i 

so that rj = \u\/V, to zeroth order in rn/M. As P{w) = 
P{'w) is isotropic we can separate out the integrals 

w = y dww^P{w) J (i^WVaoiww), (12) 

and consider different distributions of ejection speeds sep- 
arately. The integration only extends over unbound or- 
bits. For w > y^2GM/r + V, all ejecta escape, and for 
w < ^y2GM/r — V, no ejecta can escape. 

Because of reflection symmetry with respect to the bi- 
nary orbital plane, u- vanishes identically. The kick is thus 
given in the orbital plane of the pre-explosion binary. It is 
easy to show that when the particles are ejected with cqiial 
velocity Wg so that P(ui) = S{w — Wo)/{'iTrWo), the inte- 
grals yield — —{GM/r)w~^ and Uy ~ V for very large 
values of Wo (^ V) which is consistent with the require- 
ment that at large ejection speeds, the outflow reaches 
infinity with an average velocity comparable to that of 
the exploding star. Figure 1 shows as a function of 
the dimensionless explosion velocity Wg/ \J GM/r. Notice 
that there is no kick (i.e., rj vanishes) when E ^ Q (see 
eq. [5]). This occurs at a critical value of the expansion 
speed Wo/^/GM/r = \/2 - V/^/GM/r which is ~ 0.414 
for V ^ -^/GM/r to zeroth order in m/M. When E < 0, 
the mass outflow vanishes identically and the ejecta fall 
back to the remaining neutron star transferring all their 
linear momentum to it. 

4. FINAL KICK SPEED FROM THE SHORT-LIVED BINARY 
4.1. Masses and orbital separations 

The flnal kick speed imparted to the remaining neutron 

star depends crucially on the evolutionary scenario leading 
to the formation and disruption of its low mass compan- 
ion. The range of possible values of , the initial mass of 
the secondary (hereon), depends on the time the fragment 
forms and how fast its cools. Both of these complications 



distinguish the problem of "early" formation and decay of 
an evanescent proto- neutron star binary (case B) from the 
problem of "late" decay of a neutron star binary (case C) 
originally discussed by Clark & Eardley (1977), Blinnikov 
et al. (1984), and Imshennik & Popov (1998); the final 
outcome may differ from the results found by the authors 
for cold neutron stars. 

Figure 2 shows an estimate of the binary separation rorb 
at the time of formation (solid line) as a function of min, 
which we regard as an independent variable here. To ob- 
tain this estimate, we assumed that the angular momen- 
tum J^,dyn of the bar-unstable proto- neutron star with 
2rot/|W^| ~ /3dyn gocs entirely into orbital angular momen- 
tum of the binary. (Even if the lump forms in corotation, 
its angular momentum is smaller than the orbital angular 
momentum of the system, substantially so if the orbital 
separation is considerably larger than its radius.) As a 
reference value for J/3,dyn we adopt ^ 4 x lO^^gcm^s"^, 
the value given by Saijo ct al. (2000) for differentially 
rotating compact stars at ftdyn (note that this is the min- 
imum J/} since the m = 2 and m = 1 modes can grow 
at values of Trot/|H^| larger than /3dyn)- As shown in the 
Figure, lighter secondaries can be accommodated on wider 
orbits where the timescalc of orbit decay for emission of 
gravitational waves tqw = 5c^r^^^/{25QG^Mtot M rrun) is 
longer (Mtot = M + min is the total mass of the collapsing 
core). 

Along the relation rorb ^"^^in for Jorb = Jfj.dym the bold 
dots show values of mmmc at different cooling times TcooI 
according to the cooling models of Strobel et al. (1999). 
The dot-dashed lines in the sample plane correspond to 
loci of fixed values of the gravitational radiation decay 
timescale for the orbit, tqw (from the top, 1, 30, 100 and 
300 seconds, as labeled in the Figure). 

Notice that gravitational radiation timescales are gener- 
ally longer than cooling timescales, with the only exception 
occurring when mmmc ^ O.7M0. Coalescence of the two 
stars is likely to occur in this case with some mass loss that 
can only be estimated via hydrodynamical simulations. In 
the opposite case when min and rorb are such that tqw is 
longer than the cooling time TcooI, the binary survives for 
a while and below we study its evolution. 



4.2. Binary evolution, mass transfer, and maximum kicks 

In the case of late formation (case C), a phase of sta- 
ble mass exchange may set in, driven by gravitational ra- 
diation emission, which terminates when the secondary 
reaches a critical mass, rntid (Clark & Eardley 1977; Blin- 
nikov et al. 1984; Jaranowski & Krolak 1992; Bildsten & 
Cutler 1992; Imshennik & Popov 1998). Subsequently, the 
low mass star (which has a core-envelope structure) loses 
its extended halo evolving along a sequence of hydrostatic 
equilibria while the orbit widens (Blinnikov et al 1984); the 
star reaches the minimum stable neutron star mass, where- 
upon it explodes (Page 1982; Colpi, Shapiro & Teukolsky 
1989, 1991, 1993; Bhnnikov et al. 1990; Sumiyoshi et al 
1998). In the case of early formation of the binary (B) we 
show below that mass transfer is likely to be unstable, and 
we argue that because of this it may just drive the star to 
the point of exploding dynamically. 
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Consider a binary in which Roche lobe overflow is in 
efi'ect; then, the radius of the low mass star R(m) over- 
fills the Roche lobe and this determines the typical orbital 
separation at the onset of mass transfer 



r = 2.2R{m){M + m)^''ym 



1/3/^1/3 



(13) 



The dotted lines in Figure 2 give the tidal radius as a 
function of m for the two different mass-radius relations 
describing a cold and warm neutron star, respectively. A 
ring/ disk of material forms around the primary star at the 
moment of Roche filling, and accretion begins to the pri- 
mary. As the disk drains angular momentum from the or- 
bit there may not be sufficient time to return it to the orbit 
through disk-donor tidal torques (see for details Bildsten 
& Cutler 1992; Lubow & Shu 1975). If 1- / is the fraction 
of orbital angular momentum stored in the disk, the or- 
bital separation would vary with time (following Bildsten 
& Cutler 1992) as 



2m(/M2 - m^) 64G^mM{M + m) 



Mm{M + m) 



5r4 



(14) 



Fitting the values of f{m/M) computed by Hut & Paczyn- 
ski (1984) for the Roche problem, Bildsten & Cutler (1992) 
found 

/ ~ (5/3)(m/M)V3 _ (3/2)(m/M)2/3. (15) 

During Roche spillover m/R{mf oc (M + m)/r^ and the 
rate of change of r is 



d\nR{m) 1\ rh 
dlnm 3 J m 



(16) 



where rh is the mass transfer rate. Using this result in 
equation (14) for the orbital evolution, we have 



1 _ d\nR{m) _ 2{fM^ - w?) 
3 dlnm M(M + m) 

Mass transfer is stable only if 

1 d\nR{m) 2{fM^-m^) 



rh _ 64G-^mM{M + m) 
m 



5r4 



dlnm 



M{M + m) 



< 



(17) 
(18) 



given the fact that m/m < 0. 

In case (C), if we adopt the mass-radius relationship 
as in Jaranowski & Krolak (1992) valid for a cold star 
(and mass below IMq), R{m) = Ram!' /{m — moY where 
i?o = 7.5 km, mo = O.OgM©, b ^ 0.79, and a = 0.83, we 
find that the condition for stable mass transfer is 



2^ ^ 1 _l_ a + (a-6)(m/mo) 
3 m/mo — 1 



2(/M ^ - m?) 
M{M + m) 



<0. (19) 



In general there is only a finite interval for the mass m for 
stability whose boundaries are determined by the two root 
of the equation T — 0; m^^^ (the lower bound) and m,^^ 
(the upper bound). When m^^^ < m < m^^^ mass transfer 
is stable. Outside this interval stability is lost. There may 
also be no roots, in which case the fate of the light star 
depends on a number of details that we will discuss below. 



Adopting a constant value for /, and setting /? ~ 7, 
to simplify analysis, we note that one regime where the 
necessary inequality fails is at low masses, near m ~ toq. 



m. 



tid 



mo 



1 



3 



2/ - 1/3 



(20) 



(For / = 0.4 and (3 



0.8, m^^^/mo = 2.7 while if 



/ = 1, Titid/^^o — 1-5.) The other regime of instability is 
at masses well above mo, in which case In i?(m)/rf In m 
is relatively small, and m = m. 



tid 



with 



m. 



1 



23 



M 2 V 36 



1 

12" 



(21) 



m+JM 



(For / = 0.4, m+JM ~ 0.4 while for / = 1, m+j 
0.8.) Thus, if there is a disk that stores orbital angular 
momentum (i.e., / < 1), the range of masses for which 
there can be stable, Roche filling transfer shrinks consid- 
erably relative to the case / = 1 explored in Blinnikov et 
al. (1984) and Imshennik & Popov (1998). 

We have computed the roots of equation (19) numer- 
ically for different R{m). For the Jaranowski & Krolak 
(1992) parametrization of R{m) for cold neutron stars (C) 
we find the following intervals of ("^tid' "^tid) units of 
Mq: (0.14,1.08) or (0.14,1.33) for / = 1 and M = IAMq 
or I.7M0; for / given by equation (15), there are no 
roots for M = IAMq, and the roots are (0.29,0.50) for 
M = I.7M0; the stability interval (TO^j,m^j) depends 
sensitively on the logarithmic derivative of R{m). Note 
that the root m[T^ exists because of the special dependence 
of R{m) on the parameter mo. Numerical hydrodynamical 
models for the description of these semi-detached tight bi- 
naries are necessary to address the problem of mass trans- 
fer, and preliminary results points toward the existence of 
stable phases (Davies, private communication). Below, we 
will adopt the values of and m^^ for a 1.7M© neutron 
star. 

What is the fate of the binary then ? In case (C) if the 
two roots m^jj and exist, and if the initial mass of 
the low mass star is m^^ < min < mj^j(^ M), the binary 
may first go through a stage of evolution via gravitational 
radiation without mass exchange, then a phase of evolu- 
tion with mass exchange at a rate determined by the rate 
of angular momentum loss due to gravitational radiation, 
and finally become unstable when mjT^ is reached (Blin- 
nikov et al. 1984; Imshennik & Popov 1998). If, on the 
other hand, mmmc ^ 'min < "^tld' then the binary tightens 
as a consequence of gravitational radiation without mass 
loss, until m fills its Roche lobe, whereupon it becomes 
unstable. If we define mgtab to be the stellar mass at the 
onset of instability, then mgtab = f^ud ™- 
rristah = rriin if mmmc < rrun < m~ 



> m^^^, and 



Hid- 



What happens once m = to^t^ sets in has been described 
by Blinnikov et al. (1984). As already mentioned, they 
argue that although the equilibrium radius of a low mass 
neutron star may become larger than the Roche radius, the 
vast majority of its mass may still be enclosed within it 
(99% is contained inside the inner 38 km or so). The strip- 
ping of the stellar mass envelope occurs slowly enough that 
the companion evolves through a series of nearly equilib- 
rium states until rummc is attained, while the orbit widens. 
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More accurately, equation (14), in the absence of gravi- 
tational wave losses, leads to a change in orbital radius r 
as a function of mass loss equal to 



dlnr 



-2/ + 



2/m 



2r77/ 



d\nm M + m M{M + m) 

which can be integrated to yield when the mass is m 



(22) 



r 



in 



m 
m 



2/ 



exp 



2(1 - /)(min - m) 



(23) 



neglecting changes in M. In this approximation the or- 
bit widens by a factor (mgtab/'Timmc)^''^ as the mass 
decreases from the value Wstab toward TOmmc- 

The orbital radius at the point of explosion is somewhat 
larger than the radius of Roche spillover, and is accord- 
ingly 

rl/3^(2/-l/3) 

(24) 



2.2ii(mstab)Mi/3TO: 



stab 



' mmc — 



mi 



The orbital speed of the low mass neutron star when it 
explodes is 



V 



GM 



0.68Gi/2Mi/3m„ 



(25) 



implying a maximum kick speed 



:ick,max 



(mmmc + M) 



V 



m. 



(/-1/6) 
stab 



[i?(mstab)]i/^ 



(26) 

that scales as M due to the dependence of V on rmmc- 

This is a direct consequence of the hypothesis of Roche 
lobe mass transfer; the system has lost memory of the 
initial orbital separation. Also note that in Vkick.max we 
have mjiimc at numerator, even when Wstab = '^tid' When 

Tlstab ■ - ' 



'tid 



O.29M0, (or when rrigtab = rn^, 



0.0925M©) the maximum kick is Vkick.max ^ 936 km s 
(800kms-i), and / - 0.4. 

Concerning case (B), warmer proto-neutron stars seem 
to satisfy a smoother mass-radius relation. A fit to Stro- 
bel et al. (1999) data gives R{m) cx m^^/^ which would 
lead to a stability interval (mmmc, 'm-tid) °^^y / close to 
imity (note that there is only a single root in this case). 
Since disk-donor torquing may act on a timescale longer 
than the lifetime of the binary, it is likely that mass trans- 
fer is always unstable under these circumstances, and as 
soon as the light star overfills its Roche lobe at min, un- 
stable mass transfer begins and there might be time for a 
phase of common envelope evolution. In this context, we 
wish to argue that the star is driven almost immediately 
toward explosion. The cooling time TcooI ^ 50 sec (case 
B; see Figure 2) exceeds the typical time for mass trans- 
fer Tmass ^ mmmc/pCsound-R^ ^ 0.02 SCC (p ^ lG"g Cm"^ 

denotes the density, close to neutron drip, of the mass los- 
ing star in its envelope, i? 50 km the size of the Roche 
lobe and Csound is the sound speed in units of 10^° kms"-'^). 
In this case, the mass losing star encounters the insta- 
bility point mmmc(T') before being stabilized by cooling. 
We speculate that the core of the light star explodes on 
its own dynamical time (~ 0.001 sec) before coalescence 



of the two stars in completed over a orbital period of 
Porb = 27r(r;'5jjjjj,/GAf )'^/^ ^ 0.03 sec. Common evolution 
is known to be accompanied by ejection of part of the enve- 
lope enshrouding the system, and this mass loss may pro- 
vide an additional thrust to the merged object. Only hy- 
drodynamical simulations can give credit to this scenario, 
and we expect a continuous transition in the physics when 
moving from C (cold late type mode) to A (hot hydrody- 
namical mode) passing through B(warm early time mode). 
If we were to compute Mcick.max by assuming that the warm 
star explodes right at the time it fills the Roche lobe, we 
would obtain nominal speeds close to 10, 000 km s~^. The 
actual value of the kick velocity is computed below comb- 
ing gravitational bending with orbital phase averaging of 
the momentum impulse. 



4.3. Velocity phase-averaging and kick speeds 

Our estimate of l^ick for cases (C) and (B) using the 
expression for rj would be not complete without consider- 
ing that explosion does happens non instantaneously. The 
star explodes on a timescale comparable to the dynamical 
timescale, and this can be close to the time of revolution 
in the binary. In case (B) in particular this effect might 
be severe. 

The final kick speed will be diminished if the explosion 
is not instantaneous to a good approximation. The orbital 
frequency when the companion explodes is 



'^mmc — 



GM 



1/2 



{).:mG'-''-iir: 



m. 



(3/-1/2) 
stab 



(27) 



which ranges from Wmmc 17 s ^ for mstab = TOmmc 
to o^mmc — 360 s"'^ for mstab = ''^tid ~ O.29Af0 and 
M = 1.7Mq (case C), and 185 for case B with 
TOmmc = O.3M0. At the time of explosion, the light 
star overfills its Roche radius i?H.cxp ^ -R(m,stab)- Mat- 
ter ejected from the unstable star at a velocity Wq crosses 
it in a time i?R,exp/M'o and we can take RR,e^p/wo as an 
estimate of the duration of the explosion, Texp- 

The explosion is almost instantaneous as long as the 
dimensionless combination 



< 1. 



Wo 



(28) 



For mstab — O.O952M0 the dimensionless ratio is 
2,800kms~^/u'o ~ 0.1 for mstab — mmmc (case C), and 
9, 000/wokms^^ ^ 0.3 for mgtab = mmmc = 0.3 in case B, 
taking Wo ~ 30, 000 km s^^ as an example. 

To be more quantitative, suppose that the star loses 
mass at a rate rn{t) = mT{t) once it becomes unstable, 
with 



/ 

Jo 



dtT{t) = 1. 



(29) 



Consider circular orbits only, and assume that although 
mass loss may extend over many orbital periods, the orbit 
of the exploding star remains unaltered during the mass 
loss. (The latter approximation should be valid if mass 
loss is rapid enough, to zeroth order in the decreasing sec- 
ondary mass.) Adopt a fixed coordinate system defined 
at f = 0, when mass loss begins, so that x = r(0)/r and 



7 



y — V{0)IV. The momentum per unit ejected mass is 
now 

u — dtT{t)[ux{x COB ut + y sin uot)] + (30) 
Jg°° dtT{t) [+Uy{-~x smujt + ycosujt)] 

where u± = {ux ± my)/V2, e± = [x ± iy)/V2, and 

r^= / dtT{t)e'^K (31) 
Jo 

Here, {ux,Uy) are the same as were computed in § |^ for 
instantaneous mass loss (at t — 0). When the finite du- 
ration of the mass foss is accounted for, the efficiency of 
the recoil is diminished from 77 (as computed in § H for 
instantaneous mass loss) to riT^. For example, if mass is 
lost at a constant rate, then T{t) — TjxpB(Tcxp — i), and 
we find 

P ^ |sin(^mincrcxp/2)| ^ ^22) 

^mmc'Texp / 2 

if instead T(t) — rpT^pexp(— t/rcxp) then 



+ (ti-'mmcT'cxp)^ 

The correction due to phase averaging, T^, is only modest 
(and somewhat model dependent) as long as uj^mcTcxp <, 
1, but for long explosions, the efficiency is diminished by 
a factor of order (wmmc''oxp)~^ in general. 

Figure 3 shows the kick velocity Vkick of the neutron star 
that remains (of mass I.7M0; but for 1.4 M© scale the ve- 
locity with M~^/'^) as a function of the speed of the ejecta 
■Wo for the two cases (C) and (B). The solid curves refer 
to scenario (C), when / is given by the Roche value (eq. 
[15]) and mstab is equal to — 0.29^/© (lower curve 
before crossover occurring at Wo ~ 40, 000 km s~^)) and 
'7^mmc (upper curve before crossover). The dotted cure 
is for mstab = m-t^j ~ 0.14Mq. The dash-dotted line, for 
case (B), is computed starting the explosion during com- 
mon envelope, when TOmmc ^ 0.3Mq and when the orbit 
separation is Porb/Tcross — PorhiV/r) times less than the 
separation at the moment of filling the Roche lobe (eq. 
[13]), to mimic orbit decay during common envelops. We 
have included the phase-average orbit correction as given 
in equation (33). 

Very high kicks come from such complex hydrodynami- 
cal process that alternative models have not been able to 
produce. 

5. DISCUSSION 

In this paper, we have elaborated on suggestions that 
substantial neutron star recoil can result from the explo- 
sion of a low mass neutron star formed in the aftermath of 
rotating core collapse (Imshennik & Popov 1998). A note- 
worthy feature of this process is that the final kick speed is 
determined by nuclear physics. The details about the for- 
mation mechanism and hydro-dynamical effects produce 
in reality a widespread range in velocities, as the one we 
see. 



We have extended earlier work by estimating the effect 
of the finite speed of the material ejected in the explo- 
sion, including both gravitational deflection, phase aver- 
aging, and indirectly taking into account cooling effects. 
We argue that a range of neutron star recoil speeds could 
arise from disruption scenarios with different companion 
masses. The recoil speed resulting from the dissolution of 
an evanescent binary system can be around 1000 kms"-'^, 
explaining the larger values deduced observationally. It 
is remarkable that nuclear physics implies a value in the 
observational range at all, once such a binary system is 
presumed to form, a concordance that lends some support 
to the idea. 

The kicks that result from this mechanism are confined 
to the orbital plane of the evanescent neutron star binary. 
We expect the spin angular momentum vector of the rem- 
nant neutron star to be nearly, if not perfectly, aligned 
with the orbital angular momentum of the binary, and in 
turn aligned with the angular momentum vector of the col- 
lapsing iron core. Near alignment would be consistent with 
the requirements imposed by the observation of geodetic 
precession of B1913+16, where the kick is constrained to 
lie very nearly in the plane of progenitor binary, which 
was most likely perpendicular to the spins of the spun- 
up neutron star (i.e. B1913-I-16) and its pre-explosion 
companion star (Wex, Kalogera, & Kramer 2000). In 
contrast, though. X-ray observations of the Vela pulsar 
have revealed a jet parallel to its proper motion (Pavlov 
et al. 2000; Helfand, Gotthelf, & Halpern 2001), and it 
has been argued that the proper motions of both Vela 
and the Crab pulsar are closely aligned with their spin 
axes (Lai, Chernoff, & Cordes 2001). The kick mecha- 
nism studied here would not be able to account for parallel 
spin and velocity. However, we note that the proper mo- 
tions of both Vela and the Crab correspond to transverse 
speeds of 70 — 141 kms~^ and 171 kms~^ respectively, us- 
ing reasonable estimates of the distances to the pulsars 
(Lai, Chernoff, & Cordes 2001). For the alignment to be 
real, the space velocities of these two systems must lie in 
the plane of the sky. The inferred speeds of these two pul- 
sars are then considerably smaller than the characteristic 
speed arising from explosion of a low mass, tidally dis- 
rupted companion. We propose that the formation of Vela 
and the Crab did not produce an evanescent binary, and 
therefore some other mechanism was responsible for their 
spin-aligned kicks (such as those explored by Lai, Chernoff, 
& Cordes 2001). We suggest that the larger kick compo- 
nent is due to formation and disruption of an evanescent 
binary, and is perpendicular to the spin axis, and that the 
smaller kick component is associated with other less vigor- 
ous kicks that tend to align with the rotation axis, perhaps 
because of phase averaging (e.g. Spruit & Phinney 1998; 
Lai, Chernoff, & Cordes 2001). A superposition of these 
two classes of kicks would also be consistent with the re- 
quirement of nearly but not precisely spin-perpendicular 
kicks to account for the observation of geodetic preces- 
sion in B1913-I-16 (Wex, Kalogera, & Kramer 2000). If 
this idea is correct, then one expects that lower velocity 
neutron stars should have their space velocities predom- 
inantly along their spin axes, and high velocity neutron 
stars should have space velocities predominantly perpen- 
dicular to their spins, giving origin to two independent 
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(low and high velocity) distributions. A contamination 
of low velocity stars (belonging to the low velocity tail 
of the high velocity distribution) would come from those 
explosions in the evanescent binary where light bending 
and rotational averaging have been more important. This 
could give rise to a himodal distribution of velocities as is 
inferred from the observations (Arzoumanian, Chernoff, & 
Cordes 2001). 

This scenario predicts a clear signature in the neutrino 
emission, in the aftermath of core collapse: Two neu- 
trino bursts occurring several seconds or minutes after core 
bounce, i.e., after the main neutrino burst, should signal 
(i) Kelvin contraction of the condensation gathered by the 
instability and (ii) explosion (on the dynamical time) of 
the light neutron star. Neutrino emission between the two 
bursts should also come form material that accretes onto 
the neutron star that remains. Heavy r— process elements, 
debris of the explosion of the light star, should also be 
ejected and found deep in the supernova expanding shells. 



Finally, we note that the same phenomenon could occur 
if the formation of a light neutron star accompanies rotat- 
ing core collapse to a black hole. In that case, we find that 
the recoil velocity of the black hole will be proportional 

2/3 1 

to Mgjj' , which is not the simple Mgjj scaling one would 
expect for a kick mechanism that ejects the same amount 
of momentum irrespective of whether a supernova leaves 
behind a neutron star or a black hole ^. The higher kick 
implied by this scaling could be compensated by stronger 
relativistic effects of gravitational bending and unstable 
mass transfer that can lead to a much lower escape prob- 
ability of the exploding debris. 

The authors thank Melvyn Davies, John Miller and An- 
drea Possenti for useful discussions, and the Referee for 
his/her critical comments. I.W. acknowledges support 
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Fig. 1. — The parameter 77 — Vkick/Msick.max as a function of the speed of the ejecta Wo expressed in units of ^ GM/r. 
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Fig. 2. — Solid line gives the orbital separation rorb (in km) versus Win, the initial mass of the light component, in a binary 
of total mass Mtot = 1.7Mq and orbital angular momentum J/3,dyn = 4 x 10^^ g cm ^ s~^. The dashed lines are the loci 

in the rorb, "T^in plane of constant time tqw for 1,30,100,300 seconds. The filled dots denote the various values of mmmc at 
the times = 1, 20, 100, 300 seconds, according to the cooling calculation of Strobcl ct al. (1999). The dotted lines give 
the binary separation at the time of Roche spillover for a cold star with R{m) — Rqiti'' /{m — mo)° as in Jaranowski & 
Krolak, and for a warm star with R{m) as given from the fit to Strobel's data: R{m) = i?*(m/m*)~^/^ with m* = 0.3Mq 
and ii* = 48 km. 
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Fig. 3. — Solid lines show the neutron star kick velocity Vkick as a function of the speed of the ejecta Wo, for case (C) 
when mstab = fntid (lower curve before crossover) and mstab = mmmc (upper curve before crossover). The dotted line is 
for mstab = ^tid ~ 0.14^^0 for / = 1. The dash-dotted line refers to case (B) for mmmc = O.SM©; for this case velocity 
phase-averaging is an important correction. 



